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A METHOD OF TREATING THE NON-GREY ERROR 
IN TOTAL EMITTANCE MEASUREMENTS 

by 

James B. Heaney and John H. Henninger 
Goddard Space Flight Center 


INTRODUCTION 

One problem that arises in the design of satellites is the variation in material properties 
between samples chosen for study and pieces actually used in fabrication. Because of the diffi- 
culties involved in determining material properties after assembly, engineers often rely on meas- 
urements obtained from materials prior to assembly, or on data taken from samples selected as 
being representative of the material used. This problem is especially crucial for thermal control 
surfaces and solar cell arrays where system performance is strongly dependent upon the surface 
condition as defined by its solar absorptance, a, and thermal emittance, e T . 

Surface materials used in spacecraft thermal control systems are selected on the basis of 
their having solar absorptance and thermal emittance values compatible with the mission’s tem- 
perature requirements. It has been shown that highly reflecting aluminum surfaces overcoated 
with dielectric films of varying thickness can produce a wide variety of coatings, with different 
a/e T values (References 1, 2). These dielectric films are non-absorbing in the solar wavelength 
region but rely on a combination of interference and absorption effects to give high infrared 
emittances. The a is essentially independent of thickness, while e T increases steadily with in- 
creasing thickness of the dielectric film. For a certain range of thickness, e T is strongly 
thickness-dependent, so a slight change in thickness can cause a rather large change in emittance. 
This property implies that the thickness must be carefully controlled and monitored. Surfaces of 
this type exhibit a spectrally varying infrared reflectance as a result of the combination of absorp- 
tion and interference effects. This arrangement often provides a compensating form of tempera- 
ture control as a byproduct, since emittance increases with temperature when the proper dielec- 
tric film thickness is chosen. 

The emittance is a sensitive function of the dielectric film thickness because the infrared 
reflectance is determined by the thickness -controlled depth of absorption bands. This spectral 
selectivity of the infrared reflectance causes the non-grey error in devices commonly used for 
the rapid determination of total emittance, and it has been shown (Reference 3) that this error can 
be quite large. Although it is often possible to eliminate this error through proper calibration 
with samples of known emittance, the report shows that when the emittance is thickness -dependent, 
calibration is either impossible or impractical. Therefore, the magnitude of the error must be 
determined for the type of device used. The choice of instrumentation is generally intended to 
satisfy the engineer’s need to correlate laboratory sample data with measurements performed on 
actual flight hardware. Several commercially available instruments permit the e T of a spacecraft 
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surface to be rapidly determined, after assembly, in a non-destructive manner. They all involve 
the determination of total emittance and are therefore subject to the non-grey error. 


EXPERIMENTAL TECHNIQUES 

Total emittance measurements were performed with a Gier-Dunkle Model DB-100 portable 
infrared reflectometer. This device produces a weighted value of the reflectance of an opaque 
sample, from which the emittance is determined. The measurement, performed relative to high 
and low reflectance standards, is instantaneous and independent of surface temperature between 
293K and 333K, Through the use of a selective filter, the non-grey error is minimized, but (as 
shown in the following calculations) not eliminated. The operating principles and design details 
of a prototype of this device have been presented in the literature (Reference 4). However, as this 
analysis proceeds, it will be necessary to briefly describe the instruments features for clarity 
and in support of the arguments. 


The values of relative weighted reflectance are obtained from Hohlraum measurements. 

Any surface covering the sample port is alternately exposed to omni-directional radiation emitted 
by hot and cold semicylinders acting as black body cavities. The radiant flux reaching the detec- 
tor from the exposed area is given by 


where 
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cos 0 cos 9 

h ~ 2 - dA 3 d V 

r*n 


A s = illuminated area of the sample, 

A d = effective viewing aperture area of the detector, 

$g , $d = angle between the indicated surface normals and the 
path, r SD , followed by the radiant flux, and 

l s = Hohlraum radiation reflected by the surface in the 
direction r SD . 


Since the illuminated area is viewed by the detector in a near-normal direction, all param- 
eters dealt with in the following analysis are considered to be averaged over the solid angle sub- 
tended by this area. For example, when the term reflectance, p , is used it is understood that 


[( 


cos 6 e cos 6 _ 

p ( *> 2 


P = 


n cos 0g cc 

S r SD 


cos 
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where p(0) is the more correct angular -dependent parameter. For the infrared region consid- 
ered here and for the narrow viewing angle used in this device, it is safe to say that p = p(B) and 
e = e{6) . 


When a test surface covering the viewing port of the instrument’s sensing head is irradiated 
by the rotating semicylindrical cavities, an alternating signal is produced in the detector because 
the cavities are at different temperatures. The intensity of radiant flux coming from the direction 
of the sample when irradiated by the hotter cavity is 


I s H = 6 s ( T S> T ( T s) * 3 (T s . T„) T (T g , T„) a Tj + C , (1) 

where 

6 S (T S ) = near-normal emittance for a sample at temperature T s , 

P s { T s , T h ) = near-normal reflectance for a sample at temperature T s when 
irradiated by the surrounding cavity walls at T H , 

^T s , t h ) = transmittance of a polyethelene compensating filter used to modify 
the spectral distribution of radiant energy reaching the detector, 

cr = Stefan- Boltzmann constant, and 

C = a term used to include all other sources of emitted or reflected 
energy reaching the detector. 

Similarly, when the surface is exposed to the colder cavity, 


\ = e s ( T s> r <T S > T s + ( T S' T c ) T (T s , T c ) a + C . (2) 

Equations (1) and (2) show that only the reflected energy varies with this alternate irradia- 
tion by the two semicylinders. The detector amplifying system is made to respond only to the 
alternating signal arising in the reflectance terms (Reference 4), so the fluctuating portion of the 
signal is the difference between Equations (1) and (2). That is, 


K v s = I Sh 


- Ic 


*s ( T S* t h> t ( T s* V * T £ - A* < T s’ T c) T ( T s’ T c) ctT c 


(3) 


where the voltage, V s , of the output signal is proportional (K) to the energy difference. The con- 
stant K includes the detector sensitivity, which is considered to be an invariant function of all the 
explicit and implicit variable parameters of Equation (3). Equation (3) represents the case for an 
unknown sample covering the opening of the sensing head. The instrument must be calibrated by 
establishing known voltage levels to define the range of the output signal. This is done using high 
and low reflectance standards, resulting in the following intensity equations for each case: 


K V ioo = *oo ( T s' V T < T s- T „) »■ Tj - P 100 (T s . T c ) T <T g . T c ) , (4) 

K V 0 = p 0 (T s , T„) r (T s , T„) a T< - P 0 (T s , T c ) r (T g , T c ) <r T< . (5) 


3 


Here p 100 ( T s , T H ) and p Q ( T s , T H ) are the reflectances of the high and low standards, respectively, 
and V 100 and v 0 are the corresponding voltages of the output signal. 

If the high and low reflectance standards have no spectral variation in reflectance over the 
wavelength range specified by the planckian distributions for crT 4 , crT^, and o-T* y and the surface 
properties of the materials are invariant over the temperature range of T H , T c , and T s? then the 
references are "grey" reflectors (Reference 5), That is, 

^100 (T s - T) or T 4 = P 100 (T s ) cr T 4 t 

and 

p 0 <T S . T) * T 4 = Po ( T s) T 4 • 


Equations (4) and (5) then become 

K V 100 = ^100 (T s ) [* (T s . T h ) aTj-r (T s , T c ) cr T 4 ] , 

K V o = P 0 ( T s> [ r (T g . T h ) cr T 4 - t (T s , T c ) cr T 4 ]. 
The measured reflectance is given by 

( T S’ T h* T c) y 


V - v 
v s v o 


v ioo V 0 


( 6 ) 

(7) 

( 8 ) 


Substituting Equations (3), (6) and (7) into Equation (8) gives, for the measured reflectance, 

[P S (T S , T h )t( t s , T h )°-T h 4 - p s (T s , T c )r(T s , T c )<tT 4 ] - p 0 (T s ) [r(T s , T H )crT 4 - r(T s , T c )crT c 4 ] 
WTs)[ t (Ts-Th>* t h 4 - T c )crT c 4 ] - p 0 (T s )[t (T s . T„) crT 4 - r(T s . T c ) crT 4 ] 


p s( T s- T H' T c) 

This reduces to 


[p s (T s , T H )r(T s , T H )aT 4 - p s (T s , T c )r(T s , T c )crT 4 ] _ 

^ ^ f r % /rr> xl f /-r nr> x 4 —XT' *T' \ nr 4*1 


^o(T s ) 


(9) 


PS( S ’ C) [^ioo(T s ) - P 0 (T S )] [r( T s , T H )aT 4 - r( T s , T c )<rT c 4 ] [<°ioo( T s) “ ^o( T S )] 

When /o 100 (T s ) and p 0 (T s ) are accurately known from independent measurements, the volt- 
ages V 100 and V 0 can be proportionately scaled so that /0 1 OO (T S ) - 100% and ^ 0 (T S ) “ 9%. Equa- 

tion (9) then reduces to 

[p s (T s , T H )r(T s , T H )aT 4 -Pg(T s , T c )r(T s , T c )o~T^] 


^°s^s J 7 *h* T c ) 


[r(T s ,T H )orT^ -r(T s ,T c )aT^] 


( 10 ) 


If, by chance, p s (T s , T h , T c ) is also the reflectance of a "grey" surface such that 

P 8 ( Tg’ 7) cr T 4 = P s (T s ) cr T 4 

as defined previously, then Equation (10) reduces further to 

Ps ( T S’ 7 h » T c ) = p s (T s ). 


(ID 


That is, for a grey reflector, the reflectance determined in this measurement is equal to the true 
reflectance by definition. 
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Equation (11) represents a special case. The more general expression of Equation (10) in- 
cludes the possibility of a sample whose spectral reflectance is non-grey* It would be erroneous 
to assume that Equation (11) is valid for all samples and that the measured reflectance is the true 
reflectance. Such an assumption would permit the measurement to be compromised by the so- 
called Tf non-grey error." 

To include the more general "non-grey" case, the parameters of Equation (10) are replaced 
by their wavelength-dependent equivalents. Equation (10) then becomes 


/•CD /*» 

J V T s)^(Ts)E x (Th)<*- J P x (T s ) V T s )E x (T c )d\ 

P S (T S ,T H ,T C )= 0 _ — * (12) 

I - I \(T s )E x (T c )dX 

*'o J o 

where crT 4 = E^(T)d\> and E^(T) is the planckian spectral irradiance function for a blackbody 

at temperature T. 

Equation (12) can be written in the form 


(T s . T h . T c ) - 


J *00 

1 ^( T s>\( T s)^(T H )-E^(T c )]dX 

n 


r x( T s ) [E^(T H ) -E X (T C )] dX 


(13) 


This equation shows that the measured reflectance equals the integrated true spectral reflectance 
modified by the indicated weighting function. By definition, the true weighted reflectance of a sur- 
face at temperature T s is given by 


P s (T s ) = 


r 

•'n 


P K < T s> E X (V d ^ 


f 

•'n 


\ (T s ) d\ 


(H) 


The difference between the true reflectance and the measured reflectance is the measure- 
ment error, defined as 


or 


^ (°s ( T s) (T s . I'h* t c), 


5 P* = 


f 


f 


V T s> V T s> dX p x (T s ) ^ (T s ) [E k (T n )-\ (T c )] dX 


/»00 -OD 

E x (T s ) d , X ^ (T s ) [ E x (T h ) - E^ (T c ) ] d X 

*'0 ^0 


(15) 


(16) 
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A summation technique is used to evaluate the integrals of Equation (16). Equation (13) may 
be rewritten to give 


where 


( T S> T H’ T c) P Akn ( T s) q n * 

n = 1 


T AU ( T s)[E A U CT H )-E Mn (T c )] 



£W T s) [ E A^V- W T c>] 

n -1 


Here /o A ^ n , r A ^, and E A ^ are defined over a wavelength band A\ rather than at a discrete wave- 
lengths . The summation is over the range 0 1 \ loo and consequently, if m is large, each AS n 
will be small. Accuracy dictates the choice of as , since the summation must follow the wave- 
length variation of each parameter in the integral. In a similar manner, Equation (14) may be 
written: 

D1 

(V ~ ^~' I P AKn < T S> p n , 

n = 1 

where 


The measurement error, defined by Equation (16) in terms of the summations presented 
above, is 


8 P s p A\ n < T s) Pn 


P A\n q " ’ 


or 

m 

8 p s = (Ts) (Pn " "O' (17) 

n = 1 

Both p n and ^ were calculated. Equations (16) and (17) show that the measurement error 
is contained within the difference between the two weighting functions, and that if they were equal 
there would be no error. The presence of the term r^(T s ) is justified by the fact that it reduces 
the difference between the two weighting functions. In most cases the surface whose emittance is 
to be determined is at room temperature, so that T s - 300K. The cavity temperatures are con- 
trolled at t h ~ 315K and T c ~ 305K. Because the difference between E^ (300K) and [e^ (315K) - 
E^ (305K)] is largest at the shorter wavelengths, a selective filter such as black polyethelene can 
minimize this difference. The spectral transmittance r K ( t s ) of the black polyethelene compen- 
sating filter is given in Figure 1 for the wavelength range 2.5 to 40 microns (^m). E AX (T S ) and 
[ E A\ (315K) - E ax (305K)] were determined from radiation tables (Reference 6). To simplify the 
calculations, only 20A\ n (i.e., m = 20 in Equation 17) were chosen with particular emphasis given 
to the wavelength region that is important for a 300K blackbody. The resulting values of p n , q n 
and (p n - q n ) are given in Table 1. 
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TRANSMITTANCE 



WAVELENGTH (*/m) 

Figure 1. Infrared transmittance of black 
polyethelene filter. 

Table 1 


Pn, q n and p n - q n for 1 < n < 20 


n 

A A. 

m) 

p„ 

q n 

Pn -q„ 

1 

0-2 

0.0000 

0.0000 

0.0000 

2 

2-4 

0.0021 

0.0026 

-0.0005 

3 

4-6 

0.0373 

0.0665 

-0.0292 

4 

6-8 

0.1010 

0.0942 

0.0068 

5 

8-10 

0.1332 

0.1555 

-0.0223 

6 

10-12 

0.1304 

0.1633 

-0.0329 

7 

12-14 

0.1124 

0.1067 

0.0057 

8 

14-16 

0.0915 

0.0777 

0.0138 

9 

16-18 

0.0727 

0.0731 

-0.0004 

10 

18-20 

0.0575 

0.0533 

0.0042 

11 

20-22 

0.0453 

0.0398 

0.0055 

12 

22-24 

0.0360 

0.0297 

0.0063 

13 

24-26 

0.0287 

0.0250 

0.0037 

14 

26-28 

0.0232 

0.0194 

0.0038 

15 

28-30 

0.0188 

0.0158 

0.0030 

16 

30-32 

0.0153 

0.0146 

0.0007 

17 

32-34 

0.0127 

0.0069 

0.0058 

18 

34-36 

0.0106 

0.0079 

0.0027 

19 

36-38 

0.0088 

0.0082 

0.0006 

20 

38-oo 

0.0625 

0.0398 

0.0227 

20 





L 

n a 1 


1.0000 

1.0000 

0.0000 


The wavelength dependence of (p n - q n ), 
from which the non-grey error arises, is 
plotted in Figure 2. The effect of the compen- 
sating filter in reducing the difference be- 
tween the weighting functions is shown by the 
dotted curve of Figure 2 which resulted from 
removing r A ^ (T ) from the equations used to 
calculate p n and q n . Equation (17) and Fig- 
ure 2 indicate that B p s is large for a sample 
having high reflectance at those wavelengths 
where (p n - q n ) is large. If the spectral re- 
flectance is unknown beforehand, as is gen- 
erally the case, a significant error can be 
recorded. 



Figure 2. Effect of modifying filter on the spectral 
distribution of the difference between weighting 
functions. 

Since p and q are normalized func- 
tions whose integrals are unity, the sum of 
their differences is zero, i.e., 

m 

J2 (Pn-<O= 0 - 

n = 1 

Therefore, when the reflectance of a surface 
is invariant with wavelength, and p ax (T s ) is 
a constant in Equation (17), then 8 p s = 0. This 
is true for a grey sample, and the conditions 
of Equation (11) apply. 
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RESULTS 


The equations previously derived are in terms of reflectance, both spectral, (T s ), and 
total, p (T s ). The measurement error can be expressed in terms of emittance, using Kirchhoff's 
relation. Since 

^ ( T S ) = 1 - ^ CTb)- 7 *. ^Bl- 
it follows that 

8 e s <T S ) = - ( T s> (18) 

for an opaque surface. 



TRUE EMITTANCE 


Figure 3. Effect of extreme spectral 
selectivity on the accuracy of meas- 
ured total emittance. 


Figure 3 presents the difference between the meas- 
ured emittance and true emittance for two hypothetical 
cases of extreme spectral selectivity (demonstrated by 
the two inserted curves). Different emittances were ob- 
tained by varying the wavelength at which the step oc- 
curred in each reflectance spectrum. The error was 
calculated from Equations (17) and (18) with the aid of 
Table 1. It is evident from the figure that the influence 
of the modifying filter reduces the maximum non-grey 
error to a tolerable level even for these extreme cases. 

To observe how the error appears in a more practical 
case, consider the reflectance curves shown in Figure 4. 
These reflectance curves demonstrate the thickness de- 
pendence of the infrared reflectance of vapor-deposited 
aluminum coated with various thicknesses of reactively- 


deposited silicon oxide. This material has been used as a thermal control surface on many satel- 
lites, and the results of a detailed study of its properties have been published (Reference 2). The 


data for Figure 4 were obtained from specular reflectance measurements performed on a Perkin- 
Elmer 621 spectrophotometer, and were used in conjunction with Equations (17) and (18) to calcu- 


late the magnitude of the non-grey error in the total emittance measurement. The calculated error 
is shown in Figure 5 for evaporated aluminum coated with various thicknesses of silicon oxide. 


To compare the magnitude of the calculated error with the empirical error, spectral re- 
flectance data for various samples were used with Equations (14) and (18) to obtain a "true emit- 
tance." Emittance values thus obtained are not subject to the non-grey error. Table 2 compares 
total emittance values measured directly, before and after correction for the non-grey error, with 
emittance values determined from weighted spectral reflectance measurements. The Alzak data 
shown in Table 2 were obtained from a Gier-Dunkle heated cavity reflectometer used with a Beck- 
man IR-7 spectrophotometer to provide infrared spectral reflectances of diffuse surfaces. The 
Alzak coating is the result of an anodic deposition of aluminum oxide onto an aluminum surface 
(Reference 7). It employs the principle described earlier for silicon oxide -over coated aluminum 
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WAVELENGTH (/jm) 

Figure 4. Infrared reflectance of aluminum 
coated with reactively deposited silicon oxide 
of various thicknesses. 


0.030 
0.020 
0.010 
S£ t 0.000 
0.010 
- 0.020 
0.030 

0.00 1.00 2.00 3.00 4.00 

SILICON OXIDE THICKNESS (//m) 

Figure 5. Calculated error in total emittance 
as a function of silicon oxide thickness. 
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Table 2 

Near Normal Total Emittance* at 300 K 



Oxide Thickness 
(mhi) 

Measured Directly 

Emittance Deduced 

Sample 

Uncorrected 

Corrected 

from Spectral 


Emittance 

Emittance 

Measurements 

Evaporated A1 

0.17 

0.013 

0.011 

0.013 

coated with 

0.36 

0.030 

0.028 

0.025 

reactively 

1.11 

0.245 

0.221 

0.222 

deposited 

1.48 

0.385 

0.360 

0.359 

silicon oxide 

1.94 

0.455 

0.443 

0.445 


2.32 

0.510 

0.505 

0.525 


3.06 

0.555 

0.558 

0.583 


2.5 

0.680 

0.682 

0.69 

Alzak 

5.0 

0.73 

0.75 

0.78 


*The third decimal place in the total emittance data is given only to illustrate the trend of the correction 
and is indicative of the reproducibility, but not the accuracy, of the measurements. 


to give a thickness -dependent a/e T , and was used as the thermal control surface on the Orbiting 
Astronomical Observatory (OAO). 

Table 2 shows that the difference between the spectrally determined and directly deter- 
mined total emittances is reduced when the calculated correction for the non-grey error is ap- 
plied to the total emittance data. The lack of perfect agreement indicates both the need for a 
better choice of weighting function in Equation (13) and the presence of other undefined instru- 
mental errors. However, the data do empirically verify the presence of the non-grey error in 
the total emittance measurement. 
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CONCLUSION 


These results show that it is not possible to eliminate the non-grey error using calibration 
techniques with currently available instrumentation when the dielectric film thickness varies in 
an unknown manner. Calibration would require knowing the infrared spectral reflectance of the 
measured surface, and this measurement would defeat from the start the requirement for a rapid 
determination of total emittance. In general, it is either impossible or impractical to measure 
the infrared spectral reflectance of coated satellite surfaces after assembly. Therefore, the only 
recourse is to calculate the magnitude of the non-grey error as it appears in a given type of in- 
strument for a particular class of samples in order to correct the total emittance measurements 
accordingly. 

A calibrated total emittance measuring device could still be used to monitor the uniformity 
of oxide thickness over a large area, since the emittance of the oxide coated surface is thickness 
dependent, but in this application it would compete with other available optical techniques. Of 
course, if the dielectric film thickness is initially known, the only need for an independent de- 
termination of emittance is to assure that the surface properties have not changed during assem- 
bly and storage. In this case, calibration is unnecessary. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, January 28, 1971 
124-09-26-16-51 
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TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
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